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(54) Method and receiver for obtaining channel estimated values for equalising digital signals 



(57) In the case of correction methods for digital 
transmission, a high efficiency can be achieved only 
when the correction parameters are adapted relatively 
accurately to the channel. To determine the parameters, 
a channel estimation is frequently firstly carried out with 
the aid of a training sequence. The quality of this chan- 
nel estimation influences the efficiency of the overall 
system decisively. The object of the invention is to make 
available a method and a receiver (50) with the aid of 
which the quality of channel estimated values obtained 
using correlative or maximum likelihood methods can 



be improved in a simple way. The efficiency of the opti- 
mum MMSE estimating method is thereby to be ap- 
proached, there being no need for advance knowledge 
of the channel in the form of coefficient variances, and 
no intention of averaging over a plurality of receiving 
blocks, which is frequently impossible in practice. This 
object is achieved according to the invention by deter- 
mining weighting factors for the original channel esti- 
mated values which are produced by mapping the esti- 
mated values from a "soft" characteristic whose varia- 
tion is influenced, for example, by the noise variance es- 
timated in the receiver. 
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Description 

[0001] The invention relates to a method for channel estinnation in a receiver of a digital transmission system, and 
to a receiver for carrying out the method. 

5 [0002] During digital transmission via dispersive channels, for example via a mobile radio channel or via two-wire 
lines, the transmitted signal is distorted and disturbed by noise. In the case of GSM, for example, the distortions are 
produced by multi-path propagation and characteristic pulse interference of the applied modulation method. A receiver 
must therefore be able to carry out special correction or detection algorithms for recovering the transmitted data from 
the received signal. The optimum method for correcting dispersive channels is Maximum Likelihood Sequence Esti- 

10 mation (MLSE) which can be carried out by means of the Viterbi algorithm. However, in the case of long channel pulse 
responses and/or non-binary signal alphabets, the Viterbi algorithm can no longer be implemented in practice because 
of the then very high outlay. Thus, suboptimum methods must be used, for example Reduced-State Sequence Esti- 
mation (RSSE) or Decision-Feedbacl< Sequence Estimation (DFSE). It is also conceivable to use Decision-Feedback 
Equalization (DFE) or linear correction, which although being of lower efficiency than sequence estimation methods 

15 also entail less system complexity. 

[0003] All known correction methods have filter coefficients as free parameters which are to be optimized as a function 
of the transmission scenario obtaining. Direct setting of these free parameters using an adaptive, recursive method, 
for example the Least-Mean Square (LMS) algorithm, which is described, for example, by S. Haykin in the article 
"Adaptive Filter Theory", Prentice -Hall, Upper Saddle River, New Jersey, third edition, 1996, is not possible in the case 

20 of a block-based transmission with relatively short blocks such as obtains, for example, in the case of the mobile radio 
standard GSM. The reason for this is that such methods require too long a convergence time. Instead of this, in such 
scenarios estimation of the channel pulse response is undertaken with the aid of a training sequence, it being possible 
to assume that the pulse response is approximately constant inside the training sequence. A substantially smaller 
sample size is required for this than in the case of direct setting. Starting from the estimated channel pulse response, 

25 it is then possible to calculate the required filter coefficients (for the special case of MLSE, the channel estimated values 
directly represent the desired parameters). In the case of relatively strongly time-variant channels, additional adaptive 
correction of the calculated coefficients (starting solution) along the block is advantageous. 

[0004] The efficiency of the downstream correction is strongly influenced by the quality of the channel estimated 
values. 

30 [0005] Transmission using pulse amplitude modulation (RAM) via a distorting channel, which produces intersymbol 

interference (ISI), can be modelled in a time-discrete fashion in accordance with Figure 1 . The sampled received signal 
t[\<l is produced as a noise-impaired convolution of the RAM transmitted sequence a[k] with the channel pulse response 
h[kl whose length is denoted by L, by means of: 

35 

r[k]^ J^h[K] a[k - k:]^ n[k]^ (1) 



40 n[k] representing the time-discrete interference which is assumed to be of Gaussian distribution and white (the latter 
being given before sampling in the case of the use of a whitened matched filter or a general root Nyquist filter as time- 
continuous receiver in put filter). Depending on the modulation method used, the amplitude coefficients and the channel 
pulse response are either purely real or complex. For the sake of simplicity, the discrete time k = -(L-1) is assigned to 
the start of each transmitting block. For the purpose of channel estimation, a training sequence known to the receiver 

45 and of length N is transmitted between the instants k= -(L-1) and k= N-L; consequently, the received symbols /[O], r 
[1] t[N-L] can be used for estimation and are now combined to form a column vector r. 



r=[r[0] t{1] ... r[N- if (2) 

50 

((.y^ denotes transposition in this case). The first L-1 received symbols t[-{L-1)], f[-1] are ignored, since these can 
still be influenced by undefined transmitted symbols; the knowledge of L required for this purpose at the receiving end 
can be presupposed. It holds on the basis of equ. (1) that 



55 



r=Ah + n, (3) 



2 



EP 1 128 617 A1 



where 



h = [h[0] h[1]... h[L-1]] , 



(4) 



n = [n[0]n[1] ... n[N-L]]' , 



(5) 



A = 



m 

a[l] 



a[-l] 
a[0] 



a[N - L] a[N - L- 1] 



a[-a - 2)] 
a[N - 'iL + 1] 



(6) 



[0006] Use is frequently made of correlative channel estinnation, which is mentioned in DE-C 1 96 14 543, for example, 
or of the maximum likelihood (ML)channel estimation related thereto in current mobile radio systems, for example 
GSM. In the case of ML estimation, the estimated value Ai^/.f^'' '® fon^ned in accordance with the criterion 



I 

^ML = arg "^2tx p{r\A, h) , 



(7) 



p(i\A, h) denoting the conditioned probability density of the received vector rfor a given matrix A and given vector h. 
It follows from equ. (7) in conjunction with consideration of the Gaussian distribution of the vector n that 



(8) 



(•)^ representing the Hermitian transpose, and the autocorrelation matrix of the noise being given by 



(9) 



Since 

n[k] was assumed to be white, it holds that 



^nn- ^N-L+P 



(10) 



with the noise variant 



and the (N-L+^) * (A/-L+1) unit matrix Ey^./^^.-]. The result is therefore, finally, a maximum likelihood channel estimated 
value of 
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(11) 



[0007] It is assumed for what follows that there is an optimum training sequence which ensures a best possible 

efficiency of ML channel estimation. This leads to 



H 



A A= cEi_, 



(12) 



the factor c in the case of PSK (Phase Shift Keying) signal constellations 



(|«W|=i). 



as used, for example, in the case of GSM or the development EDGE (Enhanced Data Rates for GSM Evolution) being 
given by c = A/-L+1 . This leads to 



(13) 



[0008] The error variance 



is in this case the same for all positions / E {0,1 L - 1} and is 



(14) 



[0009] The quality of the channel estimation can be raised by comparison with ML estimation when there is advance 
knowledge of the channel. Reference may be made in this connection, for example, to the article "Optimal Parametric 
Feedforward Estimation of Frequency-Selective Fading Radio Channels" IEEE Trans. On Gommun., COM-42: 

1 639-1 650, February, March, April, by S.A. Fechtel and H. Meyr. In this case, a maximum, a-posteriori estimated value 
for h can be calculated, being given by 



f^MAP =arginaxp(/i|/l,r) 



(15) 



= arg max{p(r | h, ^)p(*)} 



(16) 



with the conditioned density function p{h\A,r) of h with knowledge of A and rand of the density p(h) of h, by means 
of which the advance knowledge of the channel appears in equ. (1 6). In mobile radio, a complex Gaussian distribution 

without a mean can be assumed for h, mostly to a good approximation. This also leads to the estimated value of 
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^ftft' representing the autocorrelation matrix of h and being given by: 

5 



0,,=Elhh"\ (18) 



[0010] In the case of uncorrelated channel coefficients, it holds that 

0^^ = diag(02, af,...,of.^ ), (19) 

diag(A', y,...) denoting a diagonal nnatrix with the diagonal elements x, y, and o^being the variance of the ith channel 

coefficient, where i G {0,1 L - 1}. The assumption of lack of correlation is frequently justified, since taking account 

of correlations between different coefficients mostly results in only slight gains in the channel estimation. Given the 
occurrence of white noise and an optimum PSK training sequence, equ. (17) can be transformed, taking account of 
equ. (10), equ. (19) and equ. (13), into 



25 



'MAP 



= diag 



1 + 



(N-L + 1)<T 



- 1 + 



0 



(N - L + Do-f 



1_ 

N-L + 1 



(20) 



'ML 



[0011] This means that after ML channel estimation a weighting of the ith obtained channei estimated vaiue is to be 
30 undertaken using the factor 



35 



d\i] = 



tj} + (tI /(N - L + 1) 



(21) 



[0012] In the case of channel coefficients with a Gaussian distribution, the MAP channel estimation can also be 

interpreted as channel estimation with a minimum squared error (Minimum Mean-Squared Error, MMSE). 
[0013] The error variance of the ith coefficient is now 



45 



50 



E{\h,,AO-h\i]\'}- '^-^-^^ (22) 



cr: + 



as follows according to intermediate calculation, and is thus lower by a factor of d[i] than in the case of an ML channel 

estimation, compare equ. (14). ^ 

[0014] In addition to knowledge of the ML estimated values /Ja^lMj determining the weighting factors c/[/], which are 
used to determine the MAP estimated values /iaiapE'I also requires knowledge of the noise variance o^, which can 
be estimated using ML channel estimation via simulation of the received user signal, and likewise a variance of the 

corresponding channel coefficient which is, in general, not accessible to the receiver. The coefficient variances can 
therefore not be estimated as a rule by averaging, since in mobile radio systems frequent use is made of frequency 
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hopping, or there is discontinuous transnnission. The result of both is that different user blocks are generally subject 
to different channel statistics. 

[0015] In conventional methods, the mean-squared error of ML channel estimation is improved without advance 
knowledge of the channel by setting to zero in a hard fashion in a postprocessing unit estimated channel coefficients 
which fall below a threshold dependent on the maximum value maXj /imJ'I estimated interference variance. 

Such a method is described, for example, in the article "Ein Konzept zur Kanalschatzung, Interferenzeliminierung und 
Leistungsregelung fur CDMA-Mobilfunksysteme" ["A concept for channel estimation, interference elimination and out- 
put regulation for CDIVIA mobile radio systems"], which appeared in the ITG Specialist Conference 1 994, pages 77-88, 
by B. Steiner et al. Use was made for this purpose of a "hard" characteristic in accordance with the equation 



1 



fif^MiU]) = T l + sign| h^^[i] 

2 V 



(23) 



A 



which is represented in Figure 3. In accordance with [lacuna], it is possible herein to set C= (maXj\hi^i_[i\\^)/4, as taught 
in PCT/US9 1/09479. Weighting of the channel estimated values by means of such a "hard" characteristic does lead 
to improved channel estimated values, but the quality of the channel estimated values is impaired because in the case 
of small signal-to-noise ratios cj^/cj^ gains have to be acquired at the expense of relatively large losses in the case of 
large signal-to-noise ratios. 

[001 6] It is therefore the object of the invention to provide a method and a receiver with the aid of which it is possible 
for the efficiency of the optimum MMSE channel estimation to be approximately achieved. 

[0017] The core idea of the invention is to be seen in carrying out postprocessing, which is improved by comparison 
with the prior art, of obtained channel estimated values which are required, for example, for channel correction in the 
case of trellis-based estimation methods. 

[0018] The above-named technical problem is solved by the invention, on the one hand with the aid of the method 
steps of Claim 1 . 

[0019] Firstly, channel estimated values are obtained in the receiver for the pulse response of a transmission channel 
from a received corrupted data signal, and thereafter the noise variance is estimated. Subsequently, the obtained 
channel estimated values are weighted with the aid of factors which have a continuous functional relationship with the 
obtained channel estimated values which depends on the noise variance. In other words, for the purpose of post- 
processing the channel estimated values, the method according to the invention uses a "soft" characteristic by contrast 
with the prior art, which uses a "hard" characteristic. 

[0020] Only knowledge on the channel estimated values and the noise variance is required to determine the weighting 
factors. Advance knowledge on the channel pulse response in the form of coefficient variances which need to be known 
in the case of optimum MMSE channel estimation is not required. 
[0021] Advantageous developments are the subject-matter of the subclaims. 

[0022] The continuous functional^relationship between the weighting factors and the obtained channel estimated 
values is defined by the rule = f{h[i\;o^, d[i\ being the weighting factor for the ith channel estimated value, h[i] being 
the ith channel estimated value, and being the noise variance. 

[0023] In the case of an expedient approach to a solution, the continuously functional relationship is determined by 
the equation 



d[i] = max 



LMl -C,a;/(N-L + I) 



,c. 



J 



A A 

cl[i\ being the weighting factor for the ith channel estimated value, hj^^[i\ being the ith maximum likelihood channel 
estimated value, being the noise variance, N being the length of a training sequence, L being the length of the time- 
discrete pulse response, r being a number in the range from, for example, 0 to 2, C^ being a number in the range from, 
for example, 0.5 to 1 .5, and Cg being a number in the range from 0 to 1 . 

[0024] The quality of the channel estimated values can be further improved if at least some of the weighting factors 
have various continuous functional relationships with the corresponding obtained channel estimated values. 
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[0025] If the channel estimated values are connplex, different weighting factors can be used for the real and imaginary 
parts of each channel estimated value. 

[0026] The weighted channel estimated values are subsequently fed to a correction or detection algorithm, in order 
to obtain the user data contained in the received corrupted signal. 
5 [0027] The channel estimated values are obtained in a way known per se by means of a correlative or maximum 

likelihood estimation. 

[0028] The technical problem specified above is likewise solved by means of the features of Claim 9. 
[0029] For this purpose, the receiver has a device for estimating values for the pulse response of a transmission 
channel from a received corrupted data signal, and a device for estimating the noise variance. Furthermore, connected 
10 downstream of the estimating device is a device which weights the channel estimated values with the aid of factors 
which have a continuous functional relationship with the obtained channel estimated values. In this case, the functional 
relationship depends on the noise variance. A correction or detection device serves to correct or detect the received 
corrupted data signal by means of the weighted channel estimated values. 

[0030] The invention is explained in more detail below with the aid of an exemplary embodiment in conjunction with 
^5 the attached drawings, wherein: 

Fig. 1 shows a block diagram of a digital transmission link with a receiver, without postprocessing. 

Fig. 2 shows a block diagram of a digital transmission link with a receiver according to the invention. 

Fig. 3 shows a hard characteristic for a postprocessing device in accordance with the prior art, and 

20 Fig. 4 shows a soft characteristic for a postprocessing device according to the invention. 

[0031 ] Figure 2 shows a digital transmission link having a transmission channel 20 disturbed by noise n [k] and with 
a time-discrete pulse response h[k] and a receiver 50 according to the invention at whose input received symbols r[k] 
arrive. The input of the receiver 50 serves as input of an ML channel estimating device 10, at whose output there is 

25 present the estimated value A^a^J/] which is fed to the input of a postprocessing device 30. The input of the receiver 50 
and the output of the channel estimating device 10 are connected to the inputs of a noise variance estimating device 
15 which supplies an estimated noise variance at the output. The input signal of the receiver 50 and the output 
signal /ippL/l of the postprocessing device 30, and the output signal of the noise variance estimating device 15 are 
applied to a correction or detection device 40. 

30 [0032] The mode of operation of the receiver 50 and, in particular, the mode of operation of the postprocessing device 
30 are explained in more detail below. 

[0033] The ML channel estimated values present at the output of the ML channel estimating device are postprocessed 
in the postprocessing device 30, which has a soft characteristic. The postprocessed channel estimated values are then 
fed to the correction or detection device 40. 
35 [0034] After the ML channel estimation in the ML channel estimating device 10, a simulation of the undisturbed 
received signal in the noise variance estimating device 15 is used, for example, to determine for the noise variance 



40 



an estimated value from which it is possible to determine directly an estimated value for the error variance of the 
ML channel estimation /{N-L +1). Subsequently, estimated values are calculated for the weighting factors d[/], oc- 
curring in the MMSE solution, for the ML channel estimated values, compare equ. (21). 



45 d[i] = 



kM -c,&:/{N-L+D 

' (24) 



is selected as first approach for this, it being possible to interpret the numerator on the right-hand side of equ. (24) as 
50 estimated value for o^, and to interpret the denominator as estimated value for +o^/(N- L+ 1)from equ. (21). The 
factor r in equ. (24) is a number in the range from, for example, 0 to 2, preferably in the range from 1 .4 to 2, and Ci is 
a number in the range from, for example, 0.5 to 1 .5. Subsequently, the postprocessing device 30 outputs as the result 
of postprocessing of the ML channel estimated values the value of 

kp ['] = d[i] • hj^d^l (25) 
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[0035] However, it is also possible wfien using equation (24) to encounter negative weighting factors, and this appears 
to be unfavourable since said factors are generally positive in the case of the MMSE solution. In order to avoid negative 
weighting factors, it is proposed to set the weighting factors to a freely selectable threshold value G2 when the value 
of equ. (24) is less than Cg, being a nunnber in the range from 0 to 1 . This leads to the expression: 



d{i\ = max 



V 



h^,\i]\-C,all{N-L + l) 



(26) 



In the case of amplitude coefficients with a Gaussian distribution and interference with a Gaussian distribution it is 
possible given this selection of the weighting factors to achieve substantially smaller squared errors over a wide range 
of the signal-to-noise ratio a^/o^ than in the case of the ML channel estimation. 

[0036] The decisive principle of the invention is that the weighting factors are calculated, for example, as a function 
of the ML channel estimated values, which leads to the general representation of 



A 



A 



a2. 



(27) 



A 



the function depending both on and on the individual coefficients and having a "soft" characteristic, as 

represented in Figure 4. In general, there is a soft transition from C2 to 1 . Limiting cases of C2 = 0 and/or /(AJy^J/];-) = 
1 are pennitted in principle; even negative functional values are not excluded. 
[0037] The previously customary "hard" characteristic 



(23) 



represented in Figure 3, turns out, by contrast, to be plainly more unfavourable than suitable soft characteristics, since 
gains in conjunction with low signal-to-noise ratios cj^/o^ must be acquired at the expense of relatively high losses in 
the case of high signal-to-noise ratios. 

[0038] The method implemented in the postprocessing device 30 can be applied not only to ML channel estimations, 
but also to general correlative channel estimation, in which the estimated values 



KorV\ = + n\i] 



(28) 



result, it being possible to calculate the variances of the interference n[i\, denoted by 0;^ , directly in the case of a given 
training sequence and noise variance g^. The previous results for ML channel estimation can therefore also be applied 
here; for this purpose, it is necessary in each case to replace o^/(A/ - L + 1) by ai. . 

[0039] The method permits the postprocessing of channel estimated values which are determined for a correction 
method carried out subsequently, the quality of the original estimated values found with correlative or ML channel 
estimation being substantially improved, with the result that it is possible for the efficiency of the optimum MMSE 
channel estimation to be approximately achieved. This is attended by an improvement in the quality of the correction. 



Claims 



1 . Method for channel estimation in a receiver of a digital transmission system, having the following steps: 
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A A 

a) obtaining cliannel estinnated values (hi^i_[i\, h^Qf^i]) for the pulse response of a transmission channel fronn 
a received corrupted data signal; 

b) estinnating the noise variance (a^) in the receiver; and ^ ^ 

c) weighting the obtained channel estimated values (hj^i[i\, ^co/?M) with the aid of factors (d[/]) which have a 
continuous functional relationship with the obtained channel estimated values (A^^J/], ^coflt'D which depends 



A 



on the estimated noise variance (o^). 



n' 



Method according to Claim 1 , characterized in that the continuous functional^elationship is defined by the rule 
d[i\ = f{h[i\;a^), d[i\ being the weighting factorfor the ith channel estimated value, h[i] beingtheith channel estimated 



value, and being the estimated noise variance 



Method according to Claim 2, characterized in that the continuous functional relationship is determined by the 
equation 



d[i] = max 





r 









A A 

d[i\ being the weighting factorforthe ith channel estimated value, Ajy^J/], being the ith maximum likelihood channel 
estimated value, being the estimated or previously known noise variance, N being the length of a training 
sequence, L being the length of the time-discrete pulse response, and C2 being a selectable threshold value which 
is in the range from 0 to 1 . 

Method according to Claim 3, characterized in that r is in the range from 0 to 2 and C^ is in the range from 0.5 to 1 .5. 

A 

Method according to one of Claims 1 to 4, characterized in that at least some of the weighting factors {d\_i\) have 
various continuous functional relationships with the corresponding obtained channel estimated values (hi^i_[i\, hQQp^ 

[/])■ 

Method according to one of Claims 1 to 5, characterized in that the channel estimated values are complex and 
different weighting factors (d[i]) are used for the real and imaginary parts of a channel estimated value (Aiy^J/], 

Method according to one of Claims 1 to 6, characterized in that the weighted channel estimation values are fed 
to a correction or detection algorithm. 

A A 

Method according to one of Claims 1 to 7, characterized in that the channel estimated values (/i/tfj/], /'coflt'D 
are obtained by means of a correlative or maximum likelihood estimation. 

Receiver in a digital transmission system for carrying out the method according to one of Claims 1 to 8, comprising: 

A A 

a device (1 0) for estimating values (/?ytf/.[/], ^co/rf'D foi'the pulse response of a transmission channel (20) from 
a received corrupted data signal; 

a device (15) for estimating the noise variance (a^, 

a^device^(30), connected downstream ofthe estimating device (1 0), for weighting the channel estimated values 
(^mlM' ^cofiJ\) with the aid of factors (d[/l) which have a continuous functional relationship with the obtained 
channel estimated values which depends on the estimated noise variance (o^), and 
a device for correcting or detecting (40) the received corrupted data signal by means of the weighted channel 
estimated values. 
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